DESY09-161 
hep-th / yymmnn 
October 2009 



Initial time singularities and admissible initial 
states for a system of coupled scalar fields 

Jiirgen Baack^] 
Fakultat Physik, Technische Universitat Dortmund 
D - 44221 Dortmund, Germany 

and 

Nina KevlishvilJ^I 

Deutsches Elektronen Synchrotron DESY 

D - 22603 Hamburg, Germany, 
Andronikashvili Institute of Physics, GAS 
0177 Tbilisi, Georgia 



Abstract 



We discuss the problem of initial states for a system of coupled scalar 
fields out of equilibrium in the one-loop approximation. The fields consist 
of classical background fields, taken constant in space, and quantum fluctu- 
ations. If the initial state is the adiabatic vacuum, i.e., the ground state of 
a Fock space of particle excitations that diagonalize the mass matrix, the 
energy-momentum tensor is infinite at t — 0, its most singular part behaves 
as 1/t. When the system is coupled to gravity this presents a problem that 
we solve by a Bogoliubov transformation of the naive initial state. As a 
side result we also discuss the canonical formalism and the adiabatic particle 
number for such a system. Most of the formalism is presented for Minkowksi 
space. Embedding the system and its dynamics into a flat FRW universe is 
straightforward and we briefly address the essential modifications. 
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1 Introduction 



The question of initial states in nonequilibrium quantum field theory has 



found considerable interest recently [HI2lEll3ll5lEllZlEllSlliniinilI2], on 



various grounds. As a very practical aspect it was realized in numerical simu- 
lations of quantum fields in cosmology [321 E| [TJ] , that the energy-momentum 
tensor had initial time singularities if the initial state was taken as the naive 
Fock-space vacuum. These had to be removed when coupling the field to 
gravity, e.g. in a Friedmann universe. A more speculative aspect that has 
attracted some interest recently was the question, whether the choice of ini- 
tial state can be expected to leave an imprint in the CMB spectrum. In any 
case it is a question of principle, to what extent the choice of initial state is 
constrained in an interacting theory of particles by consistency requirements. 

As a general aspect of quantum field theory, the problem of initial con- 
ditions was realized long ago. It ultimately can traced back to the fact that 
one switches on the interaction at some time t = 0. In the case of nonequi- 
librium dynamics we have to impose initial conditions for the background 
fields. In most applications the initial state was taken to be the adiabatic 
vacuum, which would be the ground state if all background fields were held 
fixed forever. When one starts to evolve the system dynamically at least the 
second derivative of the fields will be discontinuous on account of the sec- 
ond order differential equation, which for t > determines their evolution. 
In cosmology it is already the first derivative of the scale parameter which 
for t > is determined by the Friedmann equation. Such singularities have 
been noted for the first time by Stiickelberg [15], they are discussed briefly 
in the textbook of Bogoliubov and Shirkov [IB] . The phenomenon has been 
identified as a kind of 'Casimir effect' connected to the initial time surface 
by Symanzik [17] . In the context of quantum field theory out of equilibrium 
the presence of such singularities has been noticed by various authors. 

The solution proposed by Symanzik is the introduction of surface coun- 
terterms in addition to the usual "bulk" counterterms of perturbative quan- 
tum field theory. In the context of nonequilibrium quantum field theory this 
line has been pursued in Refs. [HI fTUj . 




The introduction of initial time surfaces singles out the particular time 
at which one starts the evolution to the extent that the surface counterterms 
become part of the field theory for t > 0. The approach of modifying the 
naive initial vacuum state [T5J [21 E] seems to be more pragmatic; the idea 
is to find the minimal requirements on an initial state that could arise from 
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a previous dynamical evolution. The latter aspect is discussed in [Tj|. The 
technique used for constructing such an initial state consists of finding a 
Bogoliubov transformation of the naive adiabatic vacuum. For the single- 
channel case and scalar fluctuations this has been done in Ref. [2] and for 
fermion fluctuations in Ref. [20] . These results were used in Refs. [21| ITi] in 
formulating the renormalized equations in a flat FRW universe. 

Our approach is based on a mode function formalism that has been in- 
troduced, for coupled channels, in Ref. [22]. The formalism ensures the 
conservation of energy with one-loop or Hartree quantum backreaction and 
has been renormalized along the lines of Ref. [23]. There the initial state 
was chosen to be the adiabatic vacuum based on a Fock space of particle 
excitations that diagonalize the initial mass matrix. It is this initial state 
that we will improve here. In Ref. [22] renormalization is based on a per- 
turbative expansion close to standard perturbation theory, in the same way 
as in Ref. [23]. The same expansion was used in Ref. [2] for analyzing the 
initial time singularity. This analysis can be carried over in a straightforward 
way to the case of coupled fields. Most other analyses of the singular early 
time behavior were based on the eikonal expansion. An eikonal formalism 
for coupled systems has been formulated recently [23]. We are not aware, 
however, of an eikonal expansion for coupled-channel systems. 

The quantum expansion of the fluctuation fields is formulated in such a 
way that the canonical commutation relations hold for t = 0. It was not dis- 
cussed in Ref. [22] how they continue to hold for t > 0. Though this is to be 
expected it is not entirely obvious, and in fact leads some nontrivial relations 
for the fluctuation modes which prove to be useful for our formalism. This 
is discussed in Appendix [B] Another approach to the canonical formalism 
for coupled-channel systems was introduced in Ref. [23], and this is another 
reason for verifying that the scheme of Ref. [22] works correctly. 

Though our main subject here is the choice of the initial state, with hind- 
sight of coupling the system to gravity, we take the occasion for discussing 
the concept of adiabatic particle number within our formalism. This is sug- 
gested by the fact that we have to discuss Bogoliubov transformations for 
coupled systems anyway and that particle numbers are usually defined by 
the coefficients of these transformations. 

The paper is organized as follows: In Sec. [2] we introduce the model we 
want to consider, a system of two coupled quantum fields with masses and 
a general fourth order potential. We define the decomposition into classical 
and fluctuation fields and the evolution of the fluctuations. In Sec. E] we 
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discuss the behavior of the Green's function at early times and construct a 
Bogoliubov transformation in order to reduce the leading singular behavior 
such that the leading time derivatives become finite at the intial time. In 
Sec. H]we present the expectation value of the energy-momentum tensor in 
the Bogoliubov-transformed initial state. In Sec. [5] we discuss the concept 
of adiabatic particle number. Some more technical subjects are transferred 
to the Appendices: the Bogoliubov transformation for coupled systems in 
Appendix [A] and some aspects of the canonical formalism in Appendix [B] 



2 The model 



We consider a system of two coupled scalar quantum fields with a Lagrangian 
density of the form 



C 



1 1 



ft? 



(2.1) 



where the indices i,j take the values 1 and 2. The hybrid model of inflation 
with 



C = -3. 

2 M 



1 



1 



+ ipvX&X + ~rn 2 4> 2 + -(x 2 ~ v 2 ) + -(fx 



A 



4 



(2.2) 



is of this form with A 



ii 



0,A 



A, A 



22 



am 



m 2 and m 2 



—av . 

Also some models involving supersymmetric flat directions [281 122] are of this 
type. The generalization to a general mass matrix and a general fourth order 
potential is possible, but we do not want do overburden the formalism with 
a profusion of indices. Also, the limitation to two fields is not essential. 
We separate the fields into classical fields and fluctuations via 



4>i = fi + A ■ 

The classical Lagrangian density then retains the form 

1„ 1 



£(0) 



\j 2 2 



(2.3) 



(2.4) 



while the fluctuation Lagrangian, of second order in the fluctuations, becomes 

A,-,- 



£(2) 



+ 



(2.5) 
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This can be written as 



2 r 



£ |2) = E 



i=l 



(2.6) 



with 



Ml, 
M 2 
M 



12 

2 

22 



2Ai2<£iy32 , 

+ 3A22V2 + A12V 9 ? 



(2.7) 
(2.8) 



If the field is coupled to gravity in a flat FRW universe, the fluctuation mass 
matrix takes a similar form. After conformal rescaling of fields and momenta 
(see, e.g., Refs. [301 D3]) one just has to replace 



ml + (& - -)R 
o 



(2.9) 



Here a is the scale parameter, R the Ricci scalar, and the £j are the conformal 
couplings. 

In the following we restrict ourselves to homogeneous background fields 
<Pi(t), so the mass matrix depends on time only. In the FRW universe the time 
parameter is conformal time, and we have an additional time dependence via 
a(r) and R(r). 

We separate the fluctuation mass matrix into its initial value and a "po- 
tential" V via 



Ml(t)=Ml 3 (0) + V tJ (t). 



We diagonalize the initial mass matrix by 

The eigenvectors are chosen to be real, and normalized to unity: 



(2.10) 



(2.11) 



i=l 



(2.12) 



The latin subscripts refer to the field components, as before, and the Greek 
superscripts refer to the two independent solutions of the eigenvalue equation. 
We now define a set of mode functions f"(k, t) for homogeneous background 
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field in the following way: 

(i) their time evolution is determined by 

f?{k,t) + k 2 f«(k,t)+Ml(t)f«(k,t) = ; (2.13) 

(ii) the initial conditions are specified as 

/°(M) = fi, (2-14) 

/««(*, 0) = -**WE, (2-15) 
where we have introduced the frequencies 



n aQ (k) = ^+F. (2.16) 

The functions /"(&, t) form a set of linearly independent solutions of the 
system of mode equations. 

The fields ^(x, t) are quantum fields. For a homogeneous background we 
can expand them as 

U^t) = E / (27r)32 ^ Q M k )/f(M) + 4(-k)/r(M)] e ikx . (2.17) 

The canonical commutation relations are 

a a (k),4(k')] = (27r) 3 2fi Q0 (fc)5 Q/3 (5 3 (k- k') . (2.18) 

In the following we will need the two-point functions at the coincidence limit , 
the "fluctuation integrals" 

^(t) = < 0|^(x,t)^(x,t)|0 > 
_ r j3jL 

= E/ WaoW g(M)/r^')- (2.i9) 

Here the expectation value is taken in the vacuum state of a Fock space, 
whose quanta have the initial masses m a (0). This is the "adiabatic vacuum" , 
defined by 

a Q (k)|0>=0 Va,k. (2.20) 

Of course this is not the ground state of the system, and the creation and 
annihilation operators aj,(k) and a Q (k) do not create free particles with the 
masses m^. We discuss some aspects of the canonical formalism in Appendix 
|B| in particular we establish that the fluctuation integral as defined above is 
real and symmetric in i and j, though this is not apparent on the right hand 
side of Eq (JSHl). 
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3 The initial time singularity of the Green's 
function and the modified initial state 



The quantum backreaction of the fluctuations onto the classical fields can 
be derived using the closed-time-path formalism [2U E2]- For the quantum 
field theories that we consider here, it has been formulated in several seminal 
publications [331 EH ESI [18j [36]. We do not repeat this here. If one just 
considers the one-loop quantum backreaction the relevant equations take a 
rather intuitive form. The equations of motion for the classical fields become 

(pi + m^i + Xuipl + 

+3Aii^iJ r n + XxzipiFm + 2Ai2^2-^i2 = , (3.1) 

and an analogous equation for (p%. As will be analyzed below, the fluctuation 
integrals Tij(t) are singular at t — 0, the time where we start the evolution. 
As t \ it behaves as tint. Though this represents a mathematical sin- 
gularity, it is finite and even zero at t — 0. So it will not prevent us from 
starting a numerical simulation. The singular behavior becomes a problem 
when we couple the field to gravity. The dynamics of the FRW scale factor a 
is determined by the energy-momentum tensor, which involves second time 
derivatives of the two-point function. If one analyzes the energy-momentum 
tensor, one indeed finds, near t = 0, a time dependence of the form 1/t in 
T£. This then prevents one from starting the dynamical evolution. Of course, 
even in flat space this infinity is an undesirable and unphysical feature of the 
energy-momentum tensor. 

As the energy-momentum tensor is a rather involved expression, espe- 
cially after renormalization, we first consider the fluctuation integral and 
find a way to remove its initial singularity, such that its first and second time 
derivatives at t = become finite. This requires less algebra and, as we have 
seen previously [2] , this is sufficient for making the energy-momentum tensor 
finite near t = 0. 

The fluctuation integrals are ultraviolet divergent. The divergences can 
be analyzed [23] by expanding with respect to orders in V which is equivalent 
to expanding with respect to the couplings . This allows one to remove 
the divergent parts and the dynamics is determined by the remaining finite 
parts. A closer analysis shows that, on the level of fluctuation integrals, the 
contributions of zeroth and first order in V are ultraviolet divergent. One 
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finds [23], up to first order in V, 



a 3 



(2tt) 3 

1 

1 



2U 



qO 



'JiOJjO 



JiOJjO 



(2t^ 3 

V a p{t) - V Q/3 (0) cos [{Q a0 + ^ )t] 



(f2 a o 


+ ^o) 2 




1 








1 




+ ^o) 3 



fi a0 + O^o 
V^(0)sin[(fi aO + ^ o )t] 

V Q/3 (t) - V Q/3 (0) COS [(flaQ + Slpo)t] 

dt' V aP (f) cos [(O,o + fyjoX* - t')} 



(3.2) 



Several integrations by parts have been performed in order to separate the 
high momentum power behavior. The first term in the integrand is quadrati- 
cally divergent, and the one proportional to V(t) is logarithmically divergent. 
In the process of renormalization these terms are removed and included in the 
mass and coupling constant renormalizations. The contribution proportional 
to V(0) vanishes as V(0) = by definition, see Eq. f)2.10p . The terms pro- 
portional to V(t) and V(t) are finite at all times. The nonanalytic parts are 
those proportional to V(0) and V(0). Near t = we find that the momentum 
integrals which multiply V(0) and V(0) behave as 



d 3 k 



d 3 k 



(2tt) 3 2n a0 n po (n a0 + n po ) 3 



- sin[(fi a0 + ft^o)*] 
cos[(fi Q0 + tt/3o)i\ 



1 



16n 



■t\a[(m a + mp)t] , 
; t 2 ln[(m a + mp)t] . 



So in general the first and second derivatives of the fluctuation integrals 
would be infinite at t = 0. 

As we have mentioned previously there are two methods for getting rid 
of this singular behavior: either one introduces surface counterterms or one 
modifies the initial state. Our approach is the second one, and we have for- 
mulated this modification of the initial state as a Bogoliubov transformation. 
The singular behavior is obviously related to the large momentum behavior 
of the integrand. So the modification of the initial state will constrain only 
its ultraviolet behavior. We are still free to modify it at finite momenta, or 
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with contributions that vanish sufficiently fast at large momenta, as e.g. a 
thermal initial state. 

The Bogoliubov transformation and its consequences for the fluctuation 
integral are presented in detail in Appendix |A] The general concept im- 
plies that we replace our naive initial state, the vacuum state for quanta of 
masses m^o by a transformed vacuum state, annihilated by a superposition of 
annihilation operators a a (k) and creation operators a* a (— k). The essential 
formulae are: 

(T) the definition of the transformation 








C a %(k) - S a ?al 



-k) 



(3.3) 



(ii) the definition of a new vacuum state |6 > via 
a 7 (k)|0 >= d^C^ a 



n 



of) 



a a (k) 



n 



30 



10 >= 



(3.4) 



where p = C~ 1 S is a symmetric matrix; 

(iii) the definition of a matrix M a p which is introduced via the expectation 
value of a a (k)a£(k') in the new vacuum: 



< 6|a a (k)4(k')|0 >= (27r) 3 <5 3 (k - k , )2^ a0 O^M a/3 (A;) ; 

and 

(iv) the relation between p and M 

M - pM J p ] = I , 



(3.5) 



(3.6) 

which ensures that the commutator of o a (k) and ol(k') is canonical. All 
matrices which we have introduced here depend on k = |k|. 

As derived in Appendix [AJ the fluctuation integral, when evaluated in 
the Bogoliubov-transformed vacuum, takes the form 

i%(f) =< o|^(x,t)v i (x,t)|6 > 



d 3 k 



+f«(t)ff*(t)M a p 

+fr(t)f?(t)p aK * P px M KX 
+frwt(t)p aK *M KP 



~f?(t)f?(t)p^M a 



(3.7) 
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We now have to determine p in such a way as to cancel the initial sin- 
gularities which are contained in the integral over 

/"(*)//*(*)■ A s discussed 
below Eq. (I3.2p the dangerous contributions are those involving 
V a/ g(0) sin [(fl a o + Q/3o)t] and V a/ g(0) cos [f2 a0 + Q/3 )t}. They have to be can- 
celled by the terms proportional to p and p 2 generated by the Bogoliubov 
transformation. If one considers Eqs. (13. 6p and (13. 7p one realizes that the 
determination of p seems to be marred already by the nonlinear relation be- 
tween M and p. We have to realize, however, that there is no unique choice 
for p, anyway. All we need is a cancellation of the dangerous terms at large 
momenta. These contributions are divided, in the integrand, by combinations 
of Q a Q and Qpo which asymptotically behave as k~ A and k~ 5 , respectively. So 
these terms become small asymptotically, and to get the correct asymptotic 
behavior of the matrix elements of p we can work in the linear approxi- 
mation. In this approximation we have M a p ~ 5 aj 3, p l3n M aK ~ p 13 ", and 
M K \p aK * pP x ~ 0. Furthermore, we can approximate — exp(— iQ a0 t), 
whenever it appears multiplied by p. Corrections would be of order p * V; as 
V is of order p this would be of order p 2 . 

When rewriting the fluctuation integral in the Bogoliubov-transformed 
vacuum Tij we use the approximations we have just mentioned. We further 
use the expansion of /f(t) as it appears in Eq. (13. 2p . but we remove 
the renormalization parts. We then obtain, to first order in V and p 

~ [ ^ J V 1 f" fP \ n a/3 p -i(n a0 +U go )t , n aj3*J({l a0 +n S0 )t] 



' / v on n JioJjo 



v af3 (o)sm[(n a0 + n po )t] 



+7n TTT"^ / dt ' (*') cos [(Q a0 + fyroX* - t')} } . (3.8) 

The cancellation of the terms which would produce a singularity at t — 
then requires 

Impa " = 2 Q X Q rO in ^ (0), (3.9) 

Repa " = 2n\ (n In v^ (0) • ( 3 - 10 ) 
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If the fields were independent V a p would be diagonal and we would obtain 
(omitting the indices of the diagonal elements) 

lmp -2n m (3.11) 



Rep U V(0) 



and 



W = 8^V V(0), + S <3 - 12) 
for the separate Bogoliubov transformations of the two fields. This agrees 
in the approximation of large momenta with the results for the one-field 
case, Eqs. (51) and (52) of Ref. [2]. There it was possible to remove the 
contributions proportional to V(0) and V(0) for all momenta. Here these 
terms are cancelled at large momenta only. 

With Eqs. (I3.9P and (13.101) we have obtained a solution to our problem of 
initial singularities. We have to stress that there is an infinite manifold of such 
solutions, differing, e.g., by a different choice of initial occupation numbers at 
finite momenta. They all have to share the same large momentum behavior, 
however. 

Once we have p we now must determine M, using (13.61) . without any 
approximation, because otherwise our transformation would not be canonical. 
Though p appears nonlinearly, Eq. (13.61) is simply a system of four linear 
equations for the matrix elements of M. Finally, the fluctuation integral 
( 13.71) has to be computed using the exact numerical solutions f^(t) in all 
four terms of the integrand. 

For the numerical computations it is preferable to implement the Bogoli- 
ubov transformation in a different way, by redefining the mode functions. 
For this purpose we introduce 

f^(k,t) = J2 y ^mk,t)N a , + fr(k,t)p^N K ,] , (3.13) 

where the matrix iV satisfies iV x iV = M [§. It can be determined using the 
eigenvalues and eigenvectors of M. One easily verifies that 

___ f <j 3 k ~ 

^(«) = E/p^^*.«)JT(*.«) 



3 iV is not uniquely determined, all we need is one particular matrix that satisfies this 
relation. As M is Hermitian, so is N. 
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is identical to the previous definition, Eq. (13.71) . As ff and /"* are solutions 
of the same equation of motion, Eq. (12.1 3p . so is f]. It is sufficient, therefore, 
to determine as a solution to this equation with the initial conditions 
implied by Eqs. fl3TT3l) . fl27Hj) and (12151) . 

Having presented the technical procedure we would like to add some com- 
ments concerning the interpretation. In order to do so we need to discuss 
the problem of initial states in a more general way. 

The adiabatic vacuum has often been used as an initial state for preheat- 
ing simulations, maybe on the grounds that after inflation one ends up in a 
temperature zero state, i.e. an "empty" vacuum. Indeed if this is the case, 
and if the evolution of the classical fields is very slow ( "adiabatic" ) , this can 
be considered to be a reasonable guess for an initial state. Another choice 
that may be reasonable, e.g. after thermalization and in a period of adiabatic 
evolution, is a thermal initial state. This is of course not a pure state. For a 
thermal state the fluctuation integral would be replaced by 

/d 3 k 

iw)*^ 1 ^ 'wr(M). (3.15) 

with 

N a (k) = [exp{n a (k)/T) - I]' 1 . (3.16) 

Here we have used the fact that the mass matrix and therefore the fluctuation 
Hamiltonian is diagonal in the basis at t = 0. 

If one takes into account the real evolution of the system before t — then 
neither the adiabatic vacuum nor a thermal initial state will be appropriate. 
If the system has started, before t = 0, in a pure quantum state, it cannot 
have ended up, at t = 0, in a thermal state or in any other state described 
simply by particle numbers N a (k). The quantum system can be interpreted 
as a system of independent free particle only after "decoherence" , a concept 
that has been addressed in the present context in Refs. [37J [38] . But even if 
the system has started, at an earlier time, with a mixed state, the interac- 
tion with the background field will have created a coherence in the different 
components of such a state at t = 0, and a representation in the form (13.151) 
will not be possible. 

So, if one takes into account the evolution of the system prior to t — 
then one would have to know the entire prehistory or at least the prehistory 
of a long period in order to describe the state at t — with its full quantum 
coherence. This is of course not possible unless one knows how to start the 



12 



system at an earlier time, facing then the same problem. The best one can 
hope for is that after some time the system will not remember much of its 
initial state. This is presumably the case if the background fields produce 
large quantum fluctuations at later times. 

The purpose of the Bogoliubov transformation is different. If we know the 
evolution of the background fields near t = (and by continuity this means 
also shortly before t — 0), to the extent that we know V(0) and V(0), or, 
equivalently, 0(0) and 0(0), then we have a limited information on the initial 
state. Constant background fields at t < and the adiabatic vacuum state 
as initial state would produce a singularity of the first two derivatives of the 
Green's function. The Bogoliubov transformation removes this singularity, 
or, more precisely, it reduces it to higher orders in the derivatives. In this 
way it takes a minimal account of the fact that the system is not static before 
t = 0. As we have displayed above the transformation may be considered as a 
modification of the state or of the mode functions. The new state should not 
be considered as a vacuum state. The adiabatic vacuum state remains the 
lowest energy state for a given set of background fields. The analysis of the 
high momentum behavior of the fluctuation integral simply shows that the 
system will never arrive at this state if the background fields keep changing 
with time. In reality, of course, we would rather expect the quantum state 
of the fluctuations to be an excited one, particularly at low momentum. Our 
simple requirement of continuity for the Green's function does not give us 
any information on this excited state, except at high momenta. 

The state generated by the Bogoliubov transformation applied to the adi- 
abatic vacuum is a pure state. Therefore, it cannot be described by a mixed 
state with suitable particle numbers N a (k). If for some physical motivation 
we want to start with a thermal state or some other state specified by parti- 
cle numbers we have to combine two different concepts: a mixed state made 
up of different excited Fock-space states and a pure state that ensures the 
continuity of the Green' s function. For a thermal state the particle numbers 
decrease exponentially as A; — > oo; then the discontinuity of the Green' s 
function solely arises from the vacuum contribution. Its Bogoliubov trans- 
formation is well-defined and compulsory at high momenta only. There are 
then, among many others, two pragmatic ways of defining a thermal initial 
state: (i) One defines the thermal state using for all momenta the modified 
mode functions /f (fc, t). This is not a quite a thermal state, though, as the 
modified mode functions are not eigenfunctions of id/dt. (ii) One uses the 
integrand of Eq. (I3.15P with the original mode functions f"(k,t) for low 
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momenta only, and the integrand of Eq. (I3.14p at high momenta. 



4 The energy-momentum tensor 

The energy-momentum tensor for the fluctuations = (T fluc ) Mi , in a homo- 
geneous background field is diagonal and has identical space-space compo- 
nents. It may be specified by the two independent expectation values 

t 00 = e = ^ < V>;(x, *)^i(x, t) (4.1) 
+V^(x,t)V^(x,t) + A4(t)V>i(x,t)V>i(x,t) > 

and 

f£ = e-3p = < -^(x,t)^(x,t) (4.2) 

+V^(x, t)V^(x, t) + 2Mj(^ i (x ) t)^{x, t) > . 

e is the energy density and p the pressure. t 00 and t£ can be evaluated in the 
adiabatic vacuum and in the Bogoliubov-transformed vacuum in the same 
way as the fluctuation integrals. We just present the expectation values 
in the Bogoliubov-transformed vacuum, the one in the adiabatic vacuum is 
obtained by substituting p — > and M — > I. For t 00 we obtain 

«» = / (S E sb { Re ( A ° /f + /f + M ^ 

a,/3 

+ [M aP + M KX p^*p aX ] {ft ft + k 2 ftft + Ml ft ft) } • ( 4 - 3 ) 



For the trace we find 

d 3 k 1 

a,f3 



K = I E -ir { Re V KM ™ (r^ + " 2/r/f + 2M «wl) 

[M a/3 + M KX p^*p aX ] {-ft ft + k 2 ftft + 2My?ff) } • (4.4) 



+ 

Both expressions can alternatively be rewritten in terms of the modified mode 
functions of Eq. ( 13.13!) . 
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Using the equation of motion for the fluctuations we can write 



-ftl P * + k 2 ftft + Mlftft = -^JTJ) 



(4.5) 



(4.6) 



Therefore the trace can be expressed in terms of the fluctuation integrals 
(1X71) as 



Both the energy density and the trace contain second space and time deriva- 
tives of the two-point function and this can transform the mild singularities 
found in the fluctuation integrals Tij into infinities at t = 0. In Ref. [2J it was 
found that the energy density remains finite even in the adiabatic vacuum. 
As the second derivatives only appear in the kinetic terms which are diag- 
onal this analysis remains valid for the coupled-channel case. However, the 
trace of the energy-momentum tensor W in the adiabatic vacuum contains 
the second time derivative of the fluctuation integrals Ta and this behaves 
as d 2 (t In t) /dt 2 = 1/t as t \ 0. In the transformed fluctuation integrals Tij 
we have removed the dangerous terms, and so has a finite value at t = 0. 

When the field is coupled to gravity [39] the expressions (14. ip and (14.21) 
receive some further contributions that we do not want to discuss here 
in detail. They can be written in terms of the fluctuation integrals and 
their first derivatives. The fluctuation integrals themselves are not infinite 
at t = 0. The most singular of the additional terms are proportional to 
(£j — l/6)HdJ r ii /dt and appear both in t o and t^. They behave as hit as 
t \ 0; in the Bogoliubov-transformed initial state Ta is replaced by Ta and 
then the energy-momentum tensor remains finite at t — 0. 

5 The adiabatic particle number 

The adiabatic particle number is obtained by representing the fluctuation 
field at time t in terms of the adiabatic Fock space at time t. The fluctuation 
field is given by Eq. (12.171) . The adiabatic Fock space at time t is defined in 
terms of particle excitations which are eigenstates of the mass matrix Ai 2 j(t). 
We define the eigenvectors of the mass matrix by 




(4.7) 



a 




(5.1) 



15 



we again choose them to be real and normalized via 

E = 5a " 



(5.2) 



and define Q at = \Jk 2 + ra^. We further expand the fields with respect to 
the new basis as 

d 3 k 



^i(x,t) 



E 

a 

-E 



(2tt) 3 2^ 

d?k 
(2^)32 



a at (k)/" + aL t (-k)/° 



a ai (k)/«-aL(-k)^ 



ikx 



ikx 



(5.3) 
(5.4) 



where we have chosen the initial conditions for the modes f^(k, t) in analogy 
to Eqs. (12.141) and (I2.15p . Using the field expansion the new annihilation 
operators a at (k) can be expressed as 



a at (k) 



(Fxe 



-ikx 



(5.5) 



They annihilate the adiabatic vacuum defined at time t. The relation to the 
original operators a a (k) is obtained by inserting the field expansion (12.171) 
into Eq. Q. We find 



'a 



with 



-1 



n at ff(k,t) + itf(k,t) /, 



a 

it i 



■I ii 



(5.6) 

(5.7) 
(5.8) 



Using the relations ( 1B.8j) . ( 1B.9[) and ( IB. 101) of Appendix [B] it is straightfor- 
ward to verify that this is a Bogoliubov transform, i.e., that Eqs. ( 1A.23I) and 
(1A.26j) are satisfied. In terms of the matrices C and S the adiabatic particle 
number density is given by Q, 



n a (k,t) = — L- < 0|aL(k)a Qi (k)|0 >=£s°^ 

° 7 



(5.9) 



4 No summation over a. For simplicity of presentation we assume the adiabatic vacuum 
as initial state. 
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Inserting Eq. (15.81) we obtain 



n a (k,t) 



1 



E 



1 



2 n J n"f* , C"/ C~f* COL COL 

atJi Jj Ji Jj JitJjt ■ 



(5.10) 



2a 



■at 
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This has a simple interpretation: one decomposes the energy density with 
respect to the fluctuations ff t . Then n a (k, t) is obtained by dividing the part 
corresponding to the fluctuations by the frequency a at of these fluctua- 
tions. This result is analogous to the one-channel case. 

While the definition (I5.9P is suggestive we would like to add that this 
particle number does not imply a representation of the fluctuation integral in 
the form (13.151) in terms of the mode functions f^(k, t). Indeed, if one wants 
to use the representation ( 15. 31) for calculating the fluctuation integral one gets 
nontrivial contributions from the operators a a t(k)a a t(k'), aj ti (k)a^ t (k') and 
aj rt (k)aj (t (k / ) as well, see Eqs. ( 1A.38I) . ( 1A.39I) and ( 1A.40I) . These contributions 
to the fluctuation integral are negligible for large particle numbers only, i.e., 
if the matrix elements S"" 7 are much larger than the C ai . 

6 Summary 

We have addressed here two topics of the nonequilibrium dynamics of coupled 
fields in a one-loop approximation to quantum backreaction: the problem of 
the initial time singularity in the energy-momentum tensor and the definition 
of the adiabatic particle number for a system of coupled scalar fields. Along 
with these topics we have considered Bogoliubov transformations and some 
aspects of the canonical formalism for such coupled systems. 

Our main interest, as evident from the title, were the initial time sin- 
gularities. We have been able to define a Bogoliubov transformation of the 
initial state that removes the initial time singularities in such a way that the 
energy- momentum tensor is finite in the limit t \ 0. This is important if one 
considers the evolution of such a system of fields coupled to gravity. Clearly, 
this Bogoliubov transformation is constrained only at large momenta. So 
modifications that are subleading at high momenta are still acceptable. We 
had to include a discussion of the canonical formalism for a coupled-channel 
system, as some of the results were needed in the construction of the initial 
time singularities: we had to ensure that the fluctuation integrals are real 
and symmetric in the indices, as they should be on account of their definition. 
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Both the discussion of the canonical formalism and of Bogoliubov trans- 
formations for coupled-channel systems are at the same time the basis for 
defining the adiabatic particle number density. So we have derived an ex- 
pression for this density in terms of the coupled system mode functions. 
It is analogous to the definition in the single-channel case and has a sim- 
ple intuitive interpretation. Another formulation for the adiabatic particle 
number, based in an eikonal formalism and the evolution of Bogoliubov coef- 
ficients, has been presented recently [21]. As both formalisms are canonical, 
the results should be equivalent, though it may be difficult to verify this 
analytically. 

Acknowledgments 

One of us (N.K.) thanks the Humboldt Foundation for financial support, and 
the Deutsche Elektronensynchrotron DESY, Hamburg, for hospitality. 

A The Bogoliubov transformation for a cou- 
pled system 

We first recall some basic relations for the case of a single quantum field, see, 
e.g., Ref. [ID]. The Bogoliubov transformation rotates creation into annihila- 
tion operators and vice versa, such as to preserve the canonical commutation 
relations 

[a(k),a(k')] = 0, 
[a(k),a t (k')j = (27r) 3 2^ 3 (k-k') , (A.l) 
[a t (k),a t (k')j = 0. 

Furthermore the transformation has to be chosen in such a way that the 
vacuum retains its total momentum zero and remains isotropic. The most 
general form of such a transformation is then induced by the operator 

Q = -J [q(k)a\k)a\-k) - g*(fc)o(k)a(-k)] , (A.2) 

via 

a(k) = exp(Q)a(k) exp(Q t ) = exp(Q)a(k) exp(-Q) (A.3) 
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and 

|0 >= exp(Q)|0 > . 
Here q(k) is a general complex function of k — |k|.We have 

[a(k),Q] = q(k)a^-k), 

[a\-k),Q] = q*(k)a(k). 

We have in general 

^ -IV 



i(k) exp(-Q) = exp(-Q) 



n=l 



[a(k),Q],Q]...],Q] n 



(A.4) 

(A.5) 
(A.6) 

(A.7) 



where the n-th term in the sum contains n commutators. The even commu- 
tators (n = 21) yield 

[[[[[a(k),Q],Q}...},Q} 2l = \q(k)\ 2l a(^ , 



the odd ones (n = 21 + 1) yield 

[[[[[a(k),Q],Q]...],Q] 2l+1 = \q(k)\ 2l q(k)a^-k) . 
Writing q(k) = ^(k)e l5 ^ with real constants 7 and 5 we find 



(A.8) 



(A.9) 



a(k)exp(-Q) = exp(—Q) [cosh(7)a(k) - sinh(7)e l,5 a t (-k)] = exp(-Q)a(k) . 

(A. 10) 

With these preliminaries the generalization is straightforward. We have 
two sets of creation and annihilation operators a£(k) and a a (k), where a = 
1, 2 labels the two independent solutions f°(k). We have the field expansion 

/d 3 k 1 
T^E^ {*Mf?(k, t)e^ + at (-k)/- (k, t)e~^} 
a 

(A.H) 

and the commutation relations 



a a (k),4(k')l = (27r) 3 2fi a0 ^ 3 (k-k') . 



(A.12) 



The operator Q now takes the form 



* 2 



d 3 k 



E 



1 



(2tt) 3 2 V ^H^ 



g^(fc)at (k)at(-k) - ,f ^*(A;K(kK(-k) 



(A.13) 



19 



The normalization convention introduced by writing ^/f2 a o^/3o has the ad- 
vantage of keeping the functions q a ^(k) symmetric in the indices. Indeed 
this symmetry is the only restriction on these functions; as they are complex 
we have six free parameters, which are functions of k. The symmetry arises 
from the fact that the products a a (k)a (3 (— k) and a /3 (k)a Q (— k) are equiva- 
lent. On the one hand the operators commute, and on the other hand the 
arguments k and — k may be exchanged as the integration is symmetric in 
the sign of k and the functions q a ^ only depend on |k|. An asymmetric part 
of these functions would simply be summed and integrated away. We again 
have = —Q and the transformation matrix exp(Q) is unitary. 

What does not work here, at least not in a general parameterization q af3 , is 
the explicit evaluation of the transformation of the annihilation and creation 
operators. The matrix q a ~iqi$* which appears after every second step in the 
evaluation of the multiple commutators, is given by 

q q "U n V 2 + g 12 V 2 Irf + k 22 ! 2 ) ' ( } 

It is Hermitian, in analogy to the reality of \q\ 2 in the single-channel case. It 
is diagonal in two cases: (i) q 12 = and (ii) q 11 = q 22 = 0. It is instructive 
to evaluate the transformation of a a (k) in the two cases. In the first case we 
find 

fii(k) = coshGg 11 !)^) -smh(\q n \)e il " ie{qll) a\(-k) , (A.15) 
a 2 (k) = cosh(|g 22 |)a 2 (k) -sinh(|g 22 |)e iars(g22) 4(-k) , (A.16) 

i.e., a simple Bogoliubov transformation for each channel. In the second case 
we have 

a 1 (k) = cosh(|g 12 |)a 1 (k)-siiih(|g 12 |) e i ^ ia ) A /^4(-k),(A.l^ 

V "20 

a 2 (k) = cosh(|g 12 |)a 2 (k)-sinh(|g 12 |)e iar ^ 12 \/5^a t 1 (-k),(A.18) 

V "io 

i.e., an annihilation operator in channel 1 is mixed with a creation operator 
in channel 2. 

In the general case the matrix q a ~iq~i$* is not diagonal. Still we can sum 
up the series formally, as a series of matrix products. As the exponential 

20 



5«(k) = J2\ 7T \° a " a ^ - Sa "4(~ k )} ■ (A.19) 

p v p° 



series converges well this can be done even numerically. We write 

2/30 

In terms of the matrix q alS we then have 

° a " = E Tab [(M * r] ^ ' (A ' 20) 

where gg* is the matrix product q a ^q<^* and the power series is a series of 
powers of this matrix. Further, we have 



(2n + 1) 

Instead of writing these matrices as power series in g Q/3 we can ask for the con- 
ditions on C and S that follow from the requirement that the commutation 
rules should be conserved. From 

K(k),^(k')] =0 (A.22) 

one finds 

(C ai S pl - S a ~<C^) = , (A.23) 

7 

or 

CS T = SC T = (CS T ) T . (A.24) 
Multiplying from the left with C' 1 and from the right with (C T ) _1 one finds 

C- l S = (C^SY , (A.25) 

i.e., this a symmetric matrix. 

Considering the nonvanishing commutator we find 

(C a ~<C^* - S^S^*) = 8 aP , (A.26) 

7 

or, in matrix form, 

CC f -S& = I, (A.27) 
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the obvious generalization of 

cosh 2 ( 7 ) - sinh 2 ( 7 ) = 1 . (A.28) 

Instead of having to deal with two matrices it may be more convenient to 
deal with just one: the condition that the operators a a (k) annihilate the 
vacuum 10 > reads 



>= , (A.29) 

where we have defined the matrix 

p = C- X S . (A.30) 

From Eq. f]A.25|) we see that p is a symmetric matrix. Indeed we had found 
previously that we have six free parameters for the most general transforma- 
tion. 

To begin with we compute the expectation value of a a (k)aj l (k / ) in the 
new vacuum. As vacua are homogeneous and isotropic we can write 

< 6|a a (k)oJ(k / )|6 >= (27r) 3 (5 3 (k - \s!)2^fVL a0 tt po M a(3 (k) . (A.31) 

This definition implies that M is a Hermitian matrix. Using the commutation 
relations and (1A.29I) we have 

< 0|a a (k)a) 3 (k , )|0 > (A.32) 
= (27r) 3 2fi Q0 5 3 (k - k')5 Q/3 + < 6|aj,(k> Q (k)|6 > 

= (27r) 3 2fi Q0 5 3 (k - k')<W + ,/g*^ p <V** < 0|a A (-k')4 (-k)|0 > . 
In terms of the matrix M we find 

M af3y ^n^ = s a(3 n a0 + ^/n^fy p aK pP x *M XK , (A.33) 

or 

M - pM T p ] = I . (A.34) 

This can be solved explicitly for M. Using the symmetry of p it easy to 
verify, e.g., using the series expansion in pp^ , that 

M=(I-pp ] )- 1 . (A.35) 



a 7 (k)|0 >= 



ajk)- J^^aJ(-k) 



n 



130 
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Furthermore, using S = Cp it is easy to see, using Eq. (1A.27[) . that 

CC ] = {I-pp^y 1 = M . (A.36) 

This implies that knowing p the matrix C is not determined uniquely. If one 
uses a basis in which p is diagonal, CC' is diagonal as well, and we have two 
free phases in the matrix C. A further useful identity is 



Mp = pM T 



(A.37) 



It can easily be verified using again the expansion of M interms of pp' . 
We next evaluate the expectation values of the other products: 



n 



< 0|a a (k)a^(k')|0 >= t/J^V" < 0K(k)4(-k')|6 > 

V "«o 

(2n) 3 6 3 (k + k')2y^ni p^M aK (k) , 



(A.38) 



< 6|4 (k)4(k')|0 >= J^p aK * < 0|a K (-k)4(k')|0 > 



n 



(2n) 3 5 3 (k + k'^^&iop^M^k) 



(A.39) 



< O|at(k) a/J (k0|6 >= J%^ P aK Y X < 0|a K (-k)4(-k')|0 > 



(2n) 3 S 3 (k-k')2^/n^ p^*p^M KX (k) . 



(A.40) 



In terms of the matrices p and M the fluctuation integral, evaluated in 
the Bogoliubov-transformed vacuum takes the form 



TiAt) =< 6|^(x,t)^(x,t)|6 > 



d 3 k 



£ 

a,/3 



2^Q a oQ 
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f?(t)f?(t)p^M ( 



+f?Wf{t)M aP 



+fr(t)ft(t) P aK *M KP 



(A.41) 



This is the basis for determining p, this is discussed in Sec. [3J 
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As we have performed a canonical transformation it is to be expected that 
Tij is real and symmetric in % and j , as it holds for J 7 ^ . The sum of the first 
and fourth terms in the bracket can be shown to be real and symmetric in 
i and j using the relation (I A. 371) . the symmetry of p and the Hermiticity of 
M. The sum of the second and third terms in the bracket can be rewritten, 
using Eq. (|A.34j) and the symmetry in the summation over a and (3, as 



= <W?/f + p aX M KX p^*ftft + p aX M KX p^*ftf? ■ (A.42) 



The first term on the right hand side is the one that appears in the fluctuation 
integrals T%j . Its sum over a = f3 with prefactor 1 / Q a Q is real and symmetric 
in i and j, see Eq. (1B.8|) . The sum of the second and third terms on the 
right hand side is obviously symmetric in % and j. It can be shown to be real 
as well. 

The analogy of the various matrices we have defined here with the co- 
efficients obtained in the one-channel case of Ref. j2j is given by C ^ 
cosh7, S <-» sinh7exp(i5), p <-> tanh7 exp(i5), M <->• cosh 2 7 , Mp <-> 
sinh27exp(i<5)/2, and 2M — I cosh 27. 



B Canonical formalism at t > 

In Sec. [2] we have defined the fluctuation integral 

Fij{t) = < V»(x,t)Vv(x,t) > 

d 3 k 



The expression on the right hand side does not appear to be symmetric in 
the indices % and j, and does not appear to be real. On the other hand the 
fields i/>j(x, t) and ^(x, t) should commute with each other. Furthermore, the 
commutator between ^(x, t) and ipj(y,t) is given by the same integral with 
the only modification that a factor exp[zk(y — x)] appears in the integrand. 
As these fields commute as well for arbitrary x and y, the expression 

f?(k,t)ff{k,t) (B.l) 



^ 20,0 W 
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should be real. This is not quite obvious. 

To begin with we consider the commutators in x space. If we calculate 
the time derivative of the equal time commutator between the fields we get 

^i(x,t),Vi(y,t)] = [^(x,t),^(y,t)] + [^(x,t),^-(y,t)] , (B.2) 
and this is zero if the canonical commutation relations 

[^(x,t),^-(y,t)] =»V(x-y) (B.3) 

hold at time t. If we require that this relation continues to hold we get the 
condition 

^(x,t),^-(y,t)] = ft(x,*),4-(y,i)] + Mx,t),#j(y,t)] = . (B.4) 

The second term can be expressed, using the equation of motion 

- Aipj + M jk *p k = , (B.5) 

by the field commutators; the term vanishes if these commutators take their 
canonical form at time t. We have to require that the first term vanishes: 

[^(x,t),^-(y,i)]=0 (B.6) 

at time t. If this identity shall continue to hold we have to make sure that 

^(x,t),^-(y,t)] = [^(x,t),^-(y,i)] + [^(x,t),^-(y,t)] = . (B.7) 

Using again the equation of motion and the symmetry of A4ij(t) this can be 
verified, if the nontrivial commutations ( IB. 311 hold, whereupon the scheme 
closes. 

This is of course the standard way for proving the time independence of 
the canonical commutation relations. But this analysis in x space show us 
how to proceed in proving the identity 
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that guarantees the reality and symmetry of the fluctuation integrals. In 
order for the identity (1B.8[) to hold at all times, we have to require in addition 
that the identities 

Im £o^tf(M)/T(M) = (B.9) 

a 

and 

hold independent of time, anc? i/jat i/ie?/ /jo/g? at t = 0. There is no direct 
evidence for any of these relations; we just can prove that they continue to 
hold if they hold at one time. In the one-field case one just has to prove 
that the canonical commutator is satisfied at all times and that follows from 
the conservation of the Wronskian. Here the Wronskian of the fluctuations 
is given, with our initial conditions, by 

w(f a , f) = J2 {ft ft - ft ft) = 2id aP n a0 , (B.n) 

i 

where the summation is with respect to the lower indices, while in the com- 
mutators we need summations over the upper indices, weighted with 1/D, a0 . 
The proof that the conditions Eqs. flB.8P - flB.10l) hold independent of time 



goes through in analogy to the proof in x space given above, using this time 
the equations of motion (12.131) . The time derivative (on both sides) of Eq. 
( IB. 81) holds, if the relation ( IB. 91) holds. The time derivative of Eq. ( 1B.9j) can 
be shown to hold by using the equations of motion for the fluctuations, and 
assuming that the relations (IB. 81) and (IB. 101) hold. Using again the equations 
of motion, the time derivative of the relation (IB. 101) holds if Eq. ( 1B.9j) holds. 

We still have to consider the initial time t = 0. With the initial conditions 
Eqs. ( EU) and (12TTBD the relations dB~8|) and flETTOl) hold trivially as we have 
chosen the /"(0) = to be the real eigenvectors of the mass matrix. We 
could still multiply the two eigenvectors with two different phase factors 
exp(i<5 a ) without spoiling these conditions. Eq. (IB.9j) at t — reduces to 



/j fiof 



jO U V ) 



and this is the orthogonality relation dual to 



Su , (B.12) 



fiof', 



fo = <W • (B.13) 
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So at t = all three relations are satisfied, and then will so for t > 0. 

Aside from their importance for the formalism developed here the rela- 
tions flB.8j) - flB.10j) represent useful checks for numerical simulations, along 



with the time independence of the Wronskian; we have verified this numeri- 
cally. 
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